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Abstract 

A mechanism for creating well-collimated beams of neutral particles or atoms 
with spins is studied. The consideration is accomplished for a general realistic case, 
taking into account: (i) the finiteness of a cylindrical trap where the atoms to be 
shot out, are stored; (ii) the possibility of manipulating the trap magnetic fields in 
order to form different anisotropic field configurations; (iii) the presence of gravity 
curving atomic trajectories. 
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1 Introduction 



Charged particles and ions are known to be easily manipulated by means of electromag- 
netic fields. The situation is more complicated for neutral particles and atoms. In what 
follows, we shall often employ, for short, the term "atom" for both elementary neutral 
particles as well as for composite systems, such as atoms or molecules. Neutral atoms can 
be successfully confined in different traps (see reviews [1-3]). Recently [4], neutrons were 
confined in a magnetic trap during the time of about 800 seconds, that is, around 13 min- 
utes; this lifetime being mainly limited by the neutron beta decay into an electron, proton, 
and anti-neutrino, rather than by the trap losses. The formation of directed beams of 
neutral atoms is usually done in one of the two following ways. The first way can be called 
mechanical, since it requires the usage of mechanical collimators, selecting atoms from an 
isotropic distribution, after which the selected particles are released through tubes with 
a high pressure difference between the ends, as is done in molecular beam masers [5,6]. 
Another mechanism may be named resonant, because it is based on the action of laser 
beams on resonant atomic transitions [1-3]. 

A novel mechanism for creating directed beams of neutral atoms has been recently 
advanced [7,8]. This mechanism can be termed magnetic, since it involves magnetic fields, 
together with a special initial polarization of atomic spins. A short qualitative description 
of the suggested magnetic semiconfinement of atoms is directly connected with its name: 
This is an effect, when neutral atoms, stored in a trap, begin moving in one direction, 
forming a well-collimated beam, at the same time being confined from the opposite side. 

The ability to form directed well-collimated beams of neutral atoms is very important 
for various applications. For example, this is necessary for the functioning of atom lasers 
[2,3]. This can also be used for studying the coherence of trapped Bose- Einstein conden- 
sates [3,9], for considering the relative motion of binary mixtures [10], for measuring the 
one-particle density matrix by the scattering of fast atoms [11], and so on. 

In order that a method of creating atomic beams be practical, it is necessary to check 
its feasibility under a realistic situation. And this is the aim of the present paper to 
generalize the consideration of the magnetic mechanism for creating atomic beams [7,8] 
by including the following really existing essential factors. 

(i) Trap shape. Atoms to be shot out are, first, stored in a trap having a finite size 
and a particular shape [1-3]. The formation of an atomic beam is accomplished inside 
the trap whose shape, therefore, is of importance. 

(ii) Field anisotropy. Fields governing the dynamics of atoms can be prepared in 
different configurations [1-3,12]. It is useful to include in the consideration additional 
parameters regulating the anisotropy of field configurations, which would help in choosing 
optimal conditions for the formation of beams. 

(iii) Influence of gravity. Atoms, contrary to photons, possess mass. Hence atoms are 
subject to gravity that may strongly influence their dynamics [3,13]. Therefore one has 
to accurately take into account the presence of gravity. 
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2 Semiconfining regime of motion 



Let us consider a system of atoms, each having mass mo, magnetic moment /xq, and 
spin S. Introduce the quantum-mechanical averages for the spatial position of an atom 
r = {r"}, with a — x,y, z, and for a spin S = {-S""}. Denote the gravitational acceleration 
as g = and the force due to interatomic interactions by f = {/"}• The evolution 

equations for the averages r and S, in the semiclassical approach [14,15], can be presented 
in the form 

dt^ "mo 9r«+^ +mo' ^ ^ 

where B is a magnetic field acting on an atom, and 

^ = ^ S X B . (2) 

The total magnetic field of a trap can be written as a sum 

B = Bi + B2 (3) 

of a time-invariant trapping field Bi and an alternating field B2. As a trapping field, one 
often uses the quadrupole field 

Bi = B[ {xe^ + ysy + Xze,) , (4) 

typical of magnetic traps [1-3], where Bq, are unitary vectors, > is a field gradient 
and A is an anisotropy parameter that can be varied. The alternating field 

B2 = B2 {hxex + hyGy) , (5) 

with ha = hait) and /i^ + Zi^ = 1, is a transerse field serving for the stabilization of atomic 
confinement. 

As is clear from the formulas (4) and (5), these magnetic fields describe a generalized 
TOP (time orbiting potential) trap. The fields specific for the latter follow if in Eq. (4) 
one sets A = —2 and in Eq. (5) a rotating transverse field is assumed. Wc accept here a 
more general form of the trap fields in order to stress that the considered mechanism of 
beam creation does not compulsory require the usage of the TOP traps but can be realized 
in other traps as well. Moreover, contrary to the TOP trap, with a fixed anisotropy, there 
exists a rich variety of other traps, where the field anisotropy can be varied in a wide 
diapason, with the ratio of radial to axial frequencies reaching a factor of 100. More 
details on possible field configurations can be found in reviews [1-3]. Modelling here the 
trap anisotropy by the parameter A, wc demonstrate the main idea that the semiconfining 
regime of motion can always be optimized by adjusting the corresponding configuration 
of trap fields. 

In what follows, it is convenient to pass to the dimensionless space variable r = {x, y, z} 
with the dimensionless components 



r 



y 



rto -Tto -fto 
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measured in units of the characteristic length 



Rn 



B2 



(6) 



The latter is called the death circle, which describes the position of the oscillating zero 
magnetic field. The length (6) hmits the size of the trapable cloud whose actual radius 
is slightly smaller [1-3]. Introduce the characteristic frequency of atomic motion, cui, and 
that of spin motion, CU2, by the equalities 



mo-Ro 



0U2 



h 



And let us employ the notation 

G 



g 



2 ' 



7^ = 



mo-Ro 



(7) 



(8) 



where G means the dimensionless gravitational force and 7 implies a collision rate. Atomic 
collisions can be modelled by a random variable ^ = defined by means of the 

stochastic averages 



< Ut) > = , < Ut) Ut') > = ^DJap 5{t - f) , 



(9) 



in which is a diffusion rate. 

With this notation, the evolution equation (1) takes the form of the stochastic differ- 
ential equation 

and the spin equation (2) can be transformed as 

dS 



(10) 



dt 



fill 



where A = [Aap] is an antisymmetric matrix with the elements 

A12 ^ Xz , A23^x + hx , Asi^y + hy . 

The evolution equations can be simplified by taking account of the existence of small 
parameters 



UJ2 



«1, 



UJ 



W2 



where the notation 



«1, 

d 
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UJ2 



«1, 



UJ = max 
t 



dt 



(12) 



(13) 



is used. These inequalities are usually valid for realistic traps [1-3]. Then the evolution 
equations can be reduced to the consideration of atomic motion on the center manifold 



4 



by employing the scale separation approach [13,16], which involves a generalization of 
averaging techniques [17] to stochastic equations. 

Under conditions (12), the matrix A can be treated as a quasi-invariant [13,18] with 
respect to the lastly oscillating spin S. Then the solution to Eq. (11) can be presented 
as a sum 

3 

S(t) = E C^^•(^) exp{u2a,{t)t} (14) 

over the eigenvectors hj — hj{t) of the matrix A, where aj — aj{t) are the related 
eigenvalues, and the coefficients Cj are given by the initial conditions, so that 

C, = S(0)-b,(0), 

(^12^123 - ajAsi)e., + (A12A31 + ajA23)ey + {Af^ + a'j)e. 



b 



(Ah - |a,f )2 + {Al, + \a,\'){Al, + Al)]y^ 

ai,2 = ±V^12+^23 + ^31 , O3 = . 

The guiding center for the spatial variable r is described by Eq. (10) in which one has 
to substitute the form (14) and to average the right-hand side according to the rule 



T^oo r Jo 

In this way, one obtains the equation 



lim ^^"</(r,h(t), !',*)> dt. 



^=uUF + G), (15) 



in which 



F = C3 hm - r {h%e^ + hte^ + Xhle,) dt , 

T^oo T Jo 



_ x + hl)S^, + iy + hl)S',+XzS^ 

' [(x + /lg)2 + (t/ + /lO)2 + A2z2]V2 ' 

[x + h^)e^ + (y + hy)ey + Xze^ 



[(a; + /i,)2 + (i/ + /i,)2 + AV]V2 ' 

where = ha{0) and 5*0 = S"{0). In the integration, defining the effective force F, one 
should keep in mind the occurrence of another small parameter 



« 1 , (16) 



because of which r is to be treated as a quasi-invariant with respect to h. Taking, for 
concreteness, the rotating field 

hx{t) — cos cut , hy{t) — sina;^ (17) 
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results in the effective force 



[(1 + x)SS + yS', + XzS^]{xe, + ye, + 2X'ze,) 

2[(1 + 2X + X2 + y2 + A2z2)(l + ^2 + ^2 + A2z2)]l/2 " ^ > 

We have also considered other forms of the transverse field h, provided it satisfies condi- 
tions (12), (13), and (16), which yields shghtly different formulas for the effective force 
F. However, shght variations of F do not quahtatively change the dynamics of atoms. 
Therefore in what follows, we particularize the effective force by Eq. (18). 

At this point, it is worth stressing that the force (18), caused by magnetic fields, acts 
on atoms only inside the trap, while outside of the latter the force vanishes. Hence, to 
consider the motion of atoms in the whole space, including the state when they fly out of 
the trap, it is necessary to concretize the trap shape. For instance, a cylindrical trap of 
radius R and length L, measured in units of Rq, can be characterized by one of the shape 
factors 

E{r) = l-Q{x' + y^-R^) O {\z\ - \) , S(r) = exp (- ^!±^ - |J) , (19) 

in which ©(•) is a unit step function. Thus, the general expression of the force, acting on 
atoms inside as well as outside the trap, is Eq. (18) multiphed by a shape factor from 
Eq. (19). It could be possible to introduce more refined shape factors. However, to study 
the basic role of the trap finiteness, it is sufficient to consider the simple expressions (19). 
Generally, the shape factors influence atomic motion only in the intermediate region, at 
the instant the atoms leave the trap. In this region, the effective force (18) of trapping 
fields smoothly diminishes to zero. Because of this, an exact form of a shape factor 
does not change much in the principal behaviour of semiconfined atoms whose main 
properties remain practically not altered. In our numerical calculations, we have studied 
the difference resulting from the usage of the step-like and smooth shape factors from Eq. 
(19). For the same trap parameters R and L, the difference in the atomic trajectories is 
practically unnoticeable. A slight distinction can only be noticed in the phase portraits 
for velocities: The step-like shape factor causes the appearance of characteristic kinks 
that are rounded off by the smooth shape factor. 

To complete the formulation of the problem, we have to fix the initial spin polarization 
So = S(0). As is shown in Refs. [7,8], for creating directed motion of atoms, one has to 
fix the initial spin polarization so that 

5o^ = 0, Sl = Q, S^ = S. (20) 

Such initial conditions for spin polarization can be prepared in one of the ways well 
known in quantum mechanics [19] or by a short resonant pulse, as is discussed in Ref. 
[20]. Accepting the initial condition (20) and taking account of the shape factor (19), we 
come to the expression for the effective force 

F = ^ XSuz (xe^ + yey + 2X'^ze^^ (21) 
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where u — u{r), with 



ii(r) = "^'^-^ (22) 

^ ^ ~ [(l + 2a; + a;2 + t/2 + A2^2)(l + ^2 + ^2 + _x2^2)]i/2 ■ V ; 

As is easy to notice, the force (21) is invariant under the transformation 

A ^ -A , S^-S . (23) 

Therefore, it is always possible to choose such A and 5* that A^* > 0, which is assumed in 
what follows. 

It is convenient to pass to the dimensionless time measured in units of {y/XS uji)~^. 
The return to the dimensional time is made by means of the replacement t — > y/XS coit. 
Also, let us introduce the dimensionless gravitational force 

Then the evolution equation (15) can be written as the system of equations 



(Px 1 . (Py 1 

-— = - uzx + Ono , -rrr — - 



^ = W + , (25) 

where the function u = u(r) is defined in Eq. (22), with the shape factor (19). 

We accomplished numerical solution of Eqs. (25) for different field-anisotropy param- 
eters A and for different trap radii R and lengths L. The results are presented in Figs. 1 
to 4, where the z axis is along the trap axis and the latter is inclined by the 45 degrees 
to the horizon, so that 5z — ^Sx- Initial conditions for the spatial position of atoms are 
close to the trap center Tq = 0, and initial velocities of atoms, ro, are varied in the interval 
—0.1 < Tq < 0.1. The motion x{t) is similar to y{t), because of which we show only the 
X — z cross-section of phase portraits. The pictures do not qualitatively change when 
varying the trap parameters R and L in the interval [1, 10]. But the acceleration is better 
for longer traps, with larger L, as well as for larger A. The maximal velocity that can be 
achieved is of the order of Vmax ~ X^/L. The aspect ratio Rasp =< x > / < z > is small. 
Rasp ^ 0.01, showing that the beam is stretched in the 2;-direction about 100 times larger 
than in the x-direction, that is, the beam is well coUimated. The figures present calcula- 
tions for the step-like shape factor which, to our mind, more clearly illustrates the role of 
the trap finiteness. As wc have numerically checked, the usage of the smooth factor, with 
the same trap parameters, docs not alter atomic trajectories. The sole difference appears 
in the phase portraits showing the axial to radial velocities: In the case of the step-like 
shape factor there are kinks in the intermediate region where atoms leave the trap. These 
kinks are rounded off in the case of the smooth shape factor. 

Comparing these numerical calculations with those of Ref . [8] , where gravity was not 
taken into account, we come to the following conclusions. The regime of motion, called 
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semiconfinement, remains, when atoms are confined from one side of the trap by the 
minimal value 

p 'V ' 

but are not confined from another side, escaping with acceleration in this direction. For 
the present choice of coordinates, the escape direction is along the axis z. If the gravity 
component > 0, the positive gravitational force only helps to atoms to move preferably 
in the z-direction, so that all atoms are semiconfined and escape from the trap forming 
a narrow beam along the z-a.xis. However, if 5^ < 0, then gravity hampers to escape for 
some of the atoms. The fraction of atoms that remain confined because of gravity is of 
order ^^^/^\S^\■ The main visible influence of gravity on the ejected beam is in curving 
atomic trajectories in the same way as gravity curves the trajectory of a ball shot out of 
a cannon. 

To have a feeling of the values of the characteristic quantities, met above, let us give 
estimates for the case of parameters typical of alkah atoms in magnetic traps [1-3,9,21]. 
The characteristic frequency of atomic motion is a;i ~ 10^ — 10^ s~^, that of spin motion 
is UJ2 ~ 10^ — 10^ s~^. The frequency of the transverse rotating field is a; ~ 10^ — 10^ s"-*^. 
The collision rate is 7 ~ 10 s~^. Consequently, 

^~lo-^ -~lo-^ ^~io-% ^~io-% 

^2 

from where it follows that the inequalities (12) and (16) hold true. The characteristic 
radius (6) is i?o ~ 0.1 — 0.5 cm. The dimensionless gravity components, defined in Eq. 
(24), can always be made small. Thus, for A ~ 2, 5 ~ 1, i?o ~ 0.5 cm, and cui ~ 10^ — 10^ 
s~^, substituting in Eq. (24) the gravitational acceleration g"" ~ 10^ cm/s^, we get 
^ot ~ 10"^ — 10~^. When the trap axis is directed against the gravitational force, then 
the amount of atoms that cannot escape from the trap is between 1% to 10%. Hence, the 
semiconfining regime of motion is not much spoiled by gravity. The majority of atoms 
fly out of the trap, forming a well-coUimated beam with a squeezing factor of about 100; 
the maximal velocity of atoms, for a trap of length 1 — 10 cm can reach the order of 100 
cm/s. 

The role of random pair collisions can be considered in the same way as in Ref. [8]. 
These collisions do not essentially disturb the semiconfining regime provided that 

' < 1 , D = sup{Da} • 



This condition, accepting for the diffusion rate an estimate D ~ ksT/h, where ks is the 
Boltzmann constant and T is temperature, can be presented as 

These inequalities show that the disorganizing influence of atomic collisions is suppressed 

when increasing the anisotropy parameter A or lowering temperature. For the charac- 
teristic values A^" ~ 1, coi ~ 10^ — 10^ s~^, and 7 ~ 10 s~^, the critical temperature. 
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below which the semiconfining regime appears, is Tc ~ 10~^ — 10~^ K. This imphes that 
for successful functioning of an atom cannon, trapped atoms are to be sufficiently cooled 
down. Such a cooling can be effectively realized in modern traps [1^3]. 

3 Discussion 

In this paper, we have analysed the influence on the magnetic semiconfinement of neu- 
tral atoms of three factors, trap size, field anisotropy, and gravity, which had not been 
considered in the previous works [7,8]. The principal result is that, despite the action of 
gravity, the semiconfining regime can always be preserved by appropriately choosing trap 
parameters. Consequently, the suggested magnetic semiconfinement can really serve as 
an efficient mechanism for creating well-coUimated beams of neutral particles and atoms. 

In order to correctly describe the semiconfining regime of motion, it has been necessary 
to lift the adiabatic approximation that is usually employed for considering the motion 
of trapped atoms. As is known from the general theory [17], the adiabatic approximation 
may be reasonable for treating the motion of a multifrequency dynamical system that is 
close to a stationary state. The motion of atoms, permanently confined inside a trap, 
is exactly such a kind of a stationary regime. This is why the adiabatic approximation 
works reasonably well for constantly confined atoms. Corrections to this approximation, 
in the case of a TOP trap, can be estimated as being of the order of uji/uj ~ 10~^, 
which being expressed in percent, makes at most a few percent. Atomic micromotion 
of atoms in a TOP trap has recently been studied experimentally [22,23]. The detected 
nonadiabatic effects were found to be of the order of a few percent, in agreement with the 
above estimate. Anharmonic corrections to the harmonic approximation of the effective 
TOP potential are also not large for trapped atoms moving close to the trap center [24]. 

The situation is different for atoms accelerated out of trap in the semiconfining regime 
of motion. Such a regime is far from stationary, because of which the adiabatic approxi- 
mation is not applicable in principle [17]. Also, in this regime, atoms may move far from 
the trap center, which prohibits the usage of a harmonic approximation. This is why, we 
have considered the effective trap force (21) in its complete form, containing the factor 
(22), never invoking harmonic approximations. At the same time, taking into account the 
anharmonic factor (22) not only makes the consideration more general but also reduces 
the influence of a particular shape factor on the motion of atoms. This happens because 
the role of the shape factor becomes noticeable when atoms leave the trap. But the effec- 
tive force (21), due to the anharmonic factor (22), fastly diminishes to zero after atoms 
leave the death circle whose radius is usually essentially smaller than the trap size. 

As is explained in the last paragraph of the previous section, the semiconfining regime 
is easier realized after precooling to sufficiently low temperatures. In the region of the 
estimated temperatures, Bose-Einstein condensates can be created. Atomic beams out- 
coupled from condensates remind photon beams from lasers, because of which a device 
emitting condensed atoms can be called an atom laser. The known standard way for real- 
izing an output coupler for trapped atoms is by transferring atoms from a trapped state 
into an untrapped state by means of a monochromatic resonant rf field, as has been done 
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in experiments [25-27]. This method produces a beam of atoms merely faUing down in the 
field of gravity. Such a device can hardly be called an atom laser, since the very first con- 
dition on a laser is that its output could be pointed in an arbitrary direction. A directional 
extraction of sodium atoms from a trapped condensate was demonstrated by using two 
laser beams stimulating Raman transitions [28] . More information on atom lasers can be 
found in reviews [2,3]. The mechanism for creating atomic beams, we suggest, differs from 
the mentioned above in three main points: (i) it can be realized by using only magnetic 
fields, without involving additional blocking devices or external laser beams; (ii) the beam 
direction is prescribed by the trap axis that can be oriented arbitrarily; (iii) the charac- 
teristics of the semiconfinement can be varied by choosing appropriate magnetic fields. In 
this way, the proposed mechanism can be employed for creating highly-directional beams 
from atom lasers. 

It would, certainly, be nice if the effect of semiconfinement could be observed experi- 
mentally. Since we are not experimentalists, it would not be reasonable from our side to 
plunge into a discussion of technical details of experiments. Our aim has been to demon- 
strate the principal feasibility of realizing the semiconfining regime of motion. However, 
to motivate experimentalists, it would, probably, be useful to touch several points whose 
clarification could facilitate experimental work. 

First of all, we would like to emphasize that the described effect is rather general and 
can be realized in different traps, with different magnetic fields. A dynamic TOP-hke trap 
was considered above for concreteness. As we have checked, the semiconfining regime of 
motion exists for static traps as well, provided the required initial conditions for atomic 
spins are prepared. We could imagine at least three ways to achieve the necessary initial 
spin polarization. 

One possibility could be to orient spins in the z-direction by means of a guiding 
longitudinal field, while the transverse field is yet switched off. After this, the latter field 

should be quickly switched on. The main problem here is that the trap field has to be 
switched on faster than the Larmor period In juj^- This looks to be difficult to accomplish 
experimentally. Although, what seems difficult today may be successfully performed in 
future. 

Another possibility is to prepare, first, atoms with the required z-polarization of spins 
in a trap. For instance, this can be done by using the trap of Ref. [29], which is a 
quadrupole trap with a bias field along the 2;-axis. Then the spin-polarized atoms are to 
be quickly loaded into another trap with the desired field configuration. The feasibility of 
transferring atoms from one trap to another by means of a sudden transfer, as opposed 
to slow transfer, has been discussed in Ref. [30] . 

One more way could be as follows. Assume that a trap consists of two chambers, 
upper and lower. Let atoms be confined in the upper chamber. Then, acting on these 
atoms by rf field, one transfers them into an untrapped state [25-28] . Obtaining the spin 
polarization corresponding to this untrapped state, atoms fall down because of gravity, 
and pass to the lower chamber. The field configuration of the latter is to be such that 
the spin polarization of the falling atoms be that required for the semiconfining regime of 
motion. Thus, the lower chamber could emit a quasi- continuous beam of atoms. 

The features of the emitted atomic beam depend, of course, on many technical details 
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related to particular experiments. The major properties are prescribed by the magnetic 
field configuration arranged. The trap shape factor, as is explained above, plays less 
important role, when the trap size is larger than the radius of the death circle or the 
radius of an atomic cloud in a trap. A continuous atomic beam being ejected out of the 
trap will be influenced by the trap shape slightly stronger than a short pulse. In the latter 
case, when the group of atoms, after an initial acceleration reaches the trap boundary, 
one could switch off the trap field abruptly. Such a situation would exactly correspond to 
the usage of the step-like shape factor. 

Finally, one could imagine experiments for which the directed motion of atoms, or- 
ganized in the proposed way, would not necessarily end up with atoms leaving the trap. 
Some of such possibilities are mentioned in the Introduction. Probably, the most interest- 
ing suggestion would be to study the relative motion of one atomic component through 
another, as is discussed in Ref. [7]. 
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Figure Captions 



Fig. 1. The trajectories of atoms in the x — z cross-section at the initial stage of 
acceleration, when < t < 50, for the trap parameters R = 1, L = 1, A = 20, and the 
gravity components — 0.01, Sy = 0, — —0.01. 

Fig. 2. The velocities of atoms on the x — z plane for the same conditions as in Fig. 
1. The characteristic kinks is a result of a step- like shape. 

Fig. 3. Typical trajectories during the long time interval < t < 300 for the trap 
parameters R — 1, L — 1, A = 10, and the gravity components 5x — 0.05, Sy — 0, 5z — 
-0.05. 

Fig. 4. The phase portrait for atomic velocities corresponding to the parameters of 
Fig. 3. The kinks are again due to a step-like shape factor. The general picture for a 
smooth shape factor remains the same, except that the kinks become rounded off. 
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